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1. Introduction 

The term "Casimir scaling" was introduced by the present authors, together with 
Del Debbio and Faber in ref . [1] , as a criterion for the validity of various confinement 
mechanisms. Casimir scaling refers to the fact, demonstrated numerically in refs. [2, 
3] , that there is an intermediate distance interval (up to the onset of color screening) 
where the string tension due to static color sources in color group representation r 
is approximately proportional to the quadratic Casimir C T of the representation, i.e. 

a r oc C r (1.1) 

This representation dependence is seen quite convincingly in data for the SU(3) group 
obtained by Bali and by Deldar in refs. [3]. Beyond the Casimir scaling regime, higher 
color charges are screened by gluons, and the SU(N) asymptotic string tension de- 
pends only on the transformation properties (or "N-ality" ) of the representation. 
A correct theory of confinement should be able to explain both Casimir scaling at 
intermediate distances, where in particular there exists a non-zero adjoint string ten- 
sion, and the N-ality dependence of the string tension at large distances, where the 
adjoint string tension vanishes. 

For SU(N) gauge theories with A^ > 3 there are a number of color representa- 
tions, apart from the defining representation and its conjugate, in which color charge 
cannot be screened by gluons down to charge in a lower dimensional representation. 
These unscreenable representations each correspond to the lowest dimensional SU(N) 
representation with N-ality k, represented by a Young tableau with one column of 
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k boxes. Because color screening cannot occur, it is possible that Casimir scaling 
for these "/^-representations" holds asymptotically, and not just in an intermediate 
range. If so, and if Casimir scaling were exact, the Casimir scaling prediction is that 

a k _ k(N - k) 

V x ~ iv — i 

An alternative scaling law proposed for /c-string tensions is the "Sine Law" 

(jfe sinnk/N 
ox sin7r/A^ 

which is motivated by MQCD [4]. 

In fact these two scaling laws are quantitatively not so very different, for iV = 
4, 5, 6. Numerical studies indicate that the Sine Law fits the SU(N) string tension 
data better than the Casimir Law in D=4 dimensions, while in D=3 dimensions 
the Casimir Law is the better of the two [5,6]. Neither formula is a perfect fit in 
both D=3 and D=4 dimensions. Moreover, as recently pointed out by Auzzi and 
Konishi [7], the Sine Law cannot be universal among all confining SU(N) theories, 
since it is found that there are non-universal corrections to this law, for softly broken 
M = 2 supersymmetric gauge theories. 

In the large A" limit, however, Casimir scaling is exact in SU(N) gauge theories 
due to the well-known factorization property. In this limit there is no difference at 
all between the Sine and Casimir scaling laws for the /^-representations, for k finite 
as A" — > oo. In general, for /c-string tensions with k <C N, we have from either law 
the simple prediction that 

— = k {k<^N) (1.4) 

We will refer to this limiting behavior of Casimir scaling as "/c-scaling" . It should 
be stressed that since /c-scaling follows from large A" factorization, this scaling law is 
guaranteed to hold for large A" at intermediate distances (i.e. up to color screening), 
and perhaps also asymptotically. 

The /c-scaling law, because it may apply at asymptotic distance scales, poses 
an interesting challenge to the center vortex theory of confinement. Let us consider 
SU(N) gauge theory at a very large but finite value of N. The vortex theory leads 
naturally to the conclusion that the asymptotic string tension depends only on N- 
ality (i.e. on k, for /c-string tensions), and we have suggested in previous work with 
Faber [8] that Casimir scaling at intermediate distances is due to the finite thickness 
of vortices, /c-scaling, however, may hold at arbitrarily large quark separations, and 
if so it cannot be attributed to the structure of the vortex core. The interesting ques- 
tion is then whether /c-scaling is somehow implicit in the center vortex confinement 
mechanism. 
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(1.2) 



(1.3) 



In this article we would like to offer the following observation: The distribution 
of center vortices in an SU(N) gauge theory is very likely controlled by an effective 
Zn gauge theory, for reasons explained below. This theory is quasi-local at the 
color-screening scale. As N — > oo, the Z N theory goes over to a compact U(l) gauge 
theory, and a Z N Wilson loop of N-ality k corresponds to the Wilson loop of an 
object with k units of abelian charge. The charge- dependence of string tensions in a 
U(l) gauge theory, in D=3 dimensions, was worked out some time ago by Ambj0rn 
and Greensite [9]; the same charge dependence is found in Z^ and U(l) lattice gauge 
theory at strong couplings in any dimension greater than D = 2. The dependence 
is precisely /c-scaling, and is due to the formation of k distinct flux tubes forming 
between the k quark and antiquark. Because of center dominance, the string tension 
of /c-charged loops in the effective U(l) theory is identical to the string tension of 
/c-strings in the corresponding SU(N) theory at large N. In this rather simple way, 
/c-scaling at large N emerges quite naturally from the vortex mechanism, at least in 
D=3, and probably also in D=4 dimensions. 

The following sections will elaborate on this observation, noting in particular 
the relevance of center monopoles to the dynamics at large N. We will recall and 
extend a little-known fact, pointed out many years ago by Bachas and Dashen [10], 
that thin center vortices are stable classical solutions of the Wilson SU(N) lattice 
action for any N > 4. This fact, when extended to effective and renormalization- 
group improved actions, is potentially of great interest in the context of the vortex 
confinement mechanism. Finally, we consider vortex densities in the large N limit. 
It is shown that vortex densities leading to /c-scaling are well-behaved in this limit 
and, contrary to some recent work [11] based on a dilute gas approximation, do not 
grow with N. 

2. The Center Monopole Gas at Large N 

Our confidence that the effective vortex theory is a Zn gauge theory of some kind is 
based on extensive numerical computations of center-projected observables in SU(2) 
gauge theory [12-15]. Generalizing from our experience with SU(2), we suggest that 
there is an effective Z N theory which is obtained from the original SU(N) theory 
in the following way: Begin by fixing the SU(N) lattice theory to an adjoint gauge, 
where the gauge-fixing condition depends only on the link variables in the adjoint 
representation. Such adjoint gauges preserve a remnant Z^ symmetry; examples 
include all the (direct, indirect, . . . ) variants of maximal center and laplacian center 
gauges. Each link variable can be expressed as the product = z^V^, where 
G Zn, and z^In is the closest center element to on the group manifold. We 
can then write 




(2.1) 
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where the domain TZ on the SU(N) manifold consists of group elements V satisfying 

ReTr[V] > ReTr[zV] for all z e Z N ^ 1 (2.2) 
The effective Z^ theory is then formally defined by integrating over V, i.e. 

expf-Se/zM = f n DV5[F[U A ]]A[U A ]e- s W ( 2 . 3 ) 

where the gauge-fixing condition applies to the link configuration U a in the adjoint 
representation. It is not necessary to know S e ff explicitly in order to simulate the 
effective theory numerically. Evaluation of a Z N Wilson loop in the theory defined 
by S e ff(z) is equivalent to the evaluation of a center-projected Wilson loop in SU(N) 
lattice gauge theory, fixed to an adjoint gauge. 

There are some general arguments [16], but no guarantee, that center projection 
in adjoint gauges will locate confining center vortices. For this property to hold, the 
center projection must pass (at least) three tests, concerning 

1. Vortex-limited Wilson loops — Let W( n )(C) denote a Wilson loop in the 
defining representation of SU(N) gauge theory, evaluated in the subset of con- 
figurations satisfying the constraint that in each corresponding center projected 
configuration, the Z N loop takes on the value 

Z(C) = exp[2ivrn/iV] (2.4) 

Then it is required that asymptotically, for large loops, 

W {n) (C)/W {0) (C) -> exp^vrn/iv] (2.5) 

This condition is necessary if a projected loop actually identifies vortices pierc- 
ing the full loop. 

2. Center Dominance — The string tension of projected loops of a given 
N-ality in the adjoint gauge should equal the asymptotic string tension of the 
corresponding unprojected Wilson loops. This condition tells us that the center 
vortices located by the projection actually account for the full confining force. 

3. Vortex Removal — If all center vortices are removed from each full lattice 
configuration, then the string tension should vanish. This should really be a 
consequence of the previous conditions, and constitutes a check. In practice, 
one checks whether the modified loops 

W\C) = (Tr[V(C)}) = (Z(C)Tr[U(C)}) (2.6) 

have vanishing asymptotic string tension (where U(C), V(C), Z(C) are prod- 
ucts of the full link, V-link, and projected link variables, respectively, around 
the contour C). 
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If all three conditions hold, then S e ff(z) is the effective action for the confining 
Zn flux. The existence of such an effective action is only known to be true from 
numerical experiments in the case of SU(2) lattice gauge theory, in certain Laplacian 
center gauges. 1 Center dominance appears to work for SU(3) in the original version 
of Laplacian center gauge [18]; the direct and indirect variants have not yet been 
investigated in the SU(3) case. 2 In the absence of evidence to the contrary, we will 
assume that an effective Z^ gauge theory satisfying the three conditions above can 
be defined for any SU(N) pure gauge theory via an appropriate choice of adjoint 
gauge, as in the SU(2) case. 

The effective action S e ff is certainly non-local at the lattice scale, but it should be 
quasi-local at the color-screening scale. This means that the action can be expressed 
as a sum of loops (and loop products) on the lattice, and the coefficients multiplying 
large loops in the action must be exponentially suppressed with loop area, for loops 
whose extension exceeds the color-screening scale. If S e ff were not quasi- local, it 
would be hard to understand how large Zn Wilson loops could have an area law 
falloff, or how the theory could avoid long-range correlations. 

The excitations of a Z^ gauge theory are thin center vortices and, for N > 3, 
center monopoles. A center monopole is the point (D=3) or line (D=4) of intersec- 
tion of a number of vortices, whose flux adds up to an integer multiple of 2n [20] . For 
iV very large, the Zn gauge group approximates a U(l) group, and center monopoles 
go over to the usual abelian monopoles of compact QED. Because the gauge group 
is compact, Dirac lines/sheets are invisible, and isolated monopoles are stable sad- 
dlepoints of the effective action. Since the action is non-local at the lattice scale, the 
action of the monopole configurations is not set by the lattice cutoff, but rather by 
the (color-screening) scale at which the effective theory becomes quasi-local. 

Let us now consider an isolated monopole solution of the D=3 dimensional U(l) 
effective theory. The lattice Bianchi identities, coupled with spherical symmetry at 
scales large compared to the lattice spacing, are sufficient to tell us that the field 
strength due to a monopole at the origin approximates the continuum form far from 
the source 

F™ n (x) = ^i^ + F?r c (2-7) 

where F^ lvac is the contribution of the Dirac string. A dilute gas of n m monopoles 
(q = 1) and antimonopoles (q — — 1) is a superposition of such terms 

F,u(x) = £ g B e ^ _ ° 3 + F»r\*) (2-8) 

a=l \ X X a\ 

1 In maximal center gauges, center dominance holds good only to an accuracy of about 20- 
30% [17]. In variants of Laplacian center gauges, the accuracy is better than 10% [15]. 

2 There are also some SU(3) results in the direct and indirect maximal center gauge, fixed by 
over-relaxation [19]. 
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The field strength is expressed in terms of variables over a compact range 



F^x) = A u (x + p)~ A v (x) - A^x + v) + An(x) - n < A a < n 

= + 27m where - vr < F^ < n (2.9) 

Since S e ff is compact and U(l) invariant, it depends on F^ u rather than F^, and 
Dirac lines contribute nothing to the action. Moreover, if the effective U(l) theory 
is local at some scale, then a weak-field low-wavelength expansion of the U(l) action 
in powers of the field strength and its derivatives starts out with the quadratic term, 

le - 1 „ 

S eff = const. + -J2F^(x) (2.10) 

e X 

Following Polyakov [21], and substituting the superposition ( |2.8[) into the weak- 
field approximation for the effective action, leads to Coulombic interactions between 
monopoles; i.e. to a monopole Coulomb gas. It should be noted that for a theory 
with N finite, the center monopole field cannot have the form (|2.7| ) arbitrarily far 
from the source, since there is a lower limit to the non-zero flux through a 
plaquette. In a plasma, however, the field of a monopole is screened at the Debye 
length, and the monopole Coulomb gas analysis should be valid providing the U(l) 
approximation ( |2.7j ) to the Z^ monopole field holds at least up to the screening 
length. 

It was shown in ref. [9] that a large Wilson loop with k units of abelian charge, in 
a monopole Coulomb gas, has an area-law falloff with string tension equal to k times 
the string tension of a single-charged loop. By the center dominance property, the 
string tension of loops in the effective U(l) theory are the same as the asymptotic 
string tension of loops of the same U(l) charge (i.e. N-ality for finite N) in the full 
SU(N) theory. In this way we obtain the fc-scaling property in SU(N) in D = 3 
dimensions. 

3. Classical Stability of Center Vortices at N > 4 

Before turning to the density distribution of vortices in the large-N limit, we would 
like first to discuss their classical stability in the context of lattice gauge theory. 

In the continuum, instantons are stable local minima of the action of pure SU(N) 
gauge theory. It is a remarkable but little-known fact, pointed out by Bachas and 
Dashen [10] in 1982, that an analogous result holds in lattice theory: Thin center 
vortices are stable local minima of the Wilson lattice action for any SU(N) gauge 
group with N > 4. The proof is quite trivial: Any "thin" center vortex configuration 
on a classical vacuum background is gauge equivalent, in an SU(N) lattice gauge 
theory, to a lattice configuration in which each link is a center element, i.e. 

U^x) = Z^(x)I N 
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Zp(x) = exp 



N 



(n^x) = l,2,...,N-l) 



(3.1) 



It is not hard to see that any vortex-creating singular gauge transformation, operating 
on the classical vacuum, will result in a configuration which can be transformed into 
this form. If any plaquette variable is different from unity, then it has been pierced 
by a vortex. Now consider a small deformation of such thin vortex configurations 



Tt[Al(x)] « 1 



(3.2) 



The thin vortex is a classical solution if the Wilson action is a (local) minimum for 
Vfj,(x) gauge equivalent to V^x) = In- Substituting (|3.2|) into the Wilson action for 
SU(N) gauge theory 



S = ^E( 2N -nU P }-TAUl] 



(3.3) 



with j3 = 2N/g 2 , we obtain 



2N 



(3.4) 



Writing the product of VMink variables around a plaquette in the usual way, as the 
exponential of a field strength 



V p = exp[iF p ] = I N + iF p - -F 2 + ... 



and 



Z p = e~ 



D 2irin p /N 

we get, to the leading order in the field strength F p of the deformation, 







2N 1 - cos 



/27T7Z ? 



+ COS 



/ 2nn^ 



Tr[Fp] + 0(F 3 



(3.5) 
(3.6) 

(3.7) 



From this expression it is clear that the action is minimized at 

Tr[F p 2 ] = 



(3.1 



and hence that the vortex configuration ( |3.1| ) is stable, providing that at each pla- 
quette with n p > 0, the condition 



cos 



/ 2%n^ 



> o 



(3.9) 
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is satisfied. Otherwise, the vortex corresponds to a local maximum. The above 
condition for vortex stability is satisfied iff 



n„ 1 N — n n 1 . 
-T7 < 7 or — — < - 3.10 

NANA v ' 

For N = 2 and N = 3, which are the most studied cases, the stability condition 
cannot be satisfied, and center vortices are clearly unstable at the classical level. 
Beginning, however, with n p — 1, 4 at N — 5, vortex stability is obtained already 
from the classical action. 

This simple result can of course be extended beyond the simple Wilson action, 
and therein lies its physical relevance. It is obvious that thin vortices, stable or not, 
are of no real importance at weak couplings, because they are suppressed by a factor 
of order 

Vortex Area 



exp 



9 2 



(3.11) 



and do not percolate. The configurations which are of physical interest are center 
vortices having some finite thickness in physical units. In order to investigate the 
stability of vortices of thickness d (or of lesser thickness, but including quantum 
fluctuations up to scale d) starting from a lattice action at spacing a, we imagine 
following the renormalization group approach, successively applying blocking trans- 
formations of the form 



e s'M 



= J DU 5[U - F(U)]e~ s[u] (3.12) 



where U are links on the blocked lattice, and F(U) is a blocking function. The 
transformations are repeated until a lattice action with lattice spacing d is obtained. 

The Monte Carlo Renormalization Group (MCRG), and the closely related per- 
fect action approach, aim to compute effective lattice actions at large scales. It is 
assumed that only a few contours (plaquettes, 6-link loops, 8-link loops...) are im- 
portant. Given an effective action at a given length scale d, we can then ask whether 
there exist stable center vortex solutions. Let us consider, for simplicity, the class of 
improved lattice actions which consist of plaquette plus 1x2 rectangle terms, i.e. 

S = c ]T (N — ReTr[C/(P)]) + a ^ (N - ReTr[U(R)\) (3.13) 

plaq P rec R 

This class of actions has been widely discussed in the lattice literature, and includes 
the tadpole- improved action [22], the Iwasaki action [23], and two-parameter approx- 
imations [24] to the Symanzik action [25], and the DBW2 action [26]. For c , c\ > 0, 
the stability of center vortices is trivial at N > 4, and the analysis is exactly like the 
Wilson action case. However, most improved actions of this type have C\ < 0, and 
so stability must be reconsidered. 
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We will now show that the condition for the stability of of the trivial vacuum 
Ufi(x) = In in the two parameter action (|3.13|) , which is 



c + 8ci > (3.14) 



is also sufficient to guarantee that thin vortices satisfying eq. ( |3.9|) are local mimina 
of the two-parameter action, and therefore classically stable. This stability condition 
is satisfied by all of the two-parameter improved actions [22-26] mentioned above. 

It is enough to consider whether the action of a trivial vacuum or vortex configu- 
ration can be lowered in a plane by some deformation V^x) ^ In in eq. (|3.2|) . If this 
cannot be done in a plane, where the Bianchi identity can be ignored and plaquette 
variables chosen independently to minimize the action, then the action cannot be 
lowered in a volume either, since the action in a volume is just the action in a set 
of planes. Then, since we are considering the case that c% is negative, we want to 
restrict our attention to configurations in which the rectangle contributions in eq. 
Q3.13D are as large as they can possibly be, for a specified set of plaquette terms 



N — ReTr[[7p] in the plane. Consider any rectangle R containing plaquettes pi, P2, 
and write 

U(p 1 ) = Z(pi) exp[i/iei • L] 

= Z(p x )(l N + %he\L a - l -f x e\e\L a L h + . . .) 
U(p 2 ) = Z(p 2 ) exp[if 2 e 2 ■ L] 

= Z{p 2 ) (l N + if 2 e a 2 L a - l -f 2 e a 2 e b 2 L a L b + . . .) (3.15) 

where e% t2 are unit vectors, and f\, f 2 are positive. Then the rectangle term has the 
form 

U(R) = U( Pl )gU(p 2 )gi (3.16) 

where g is an SU(N) group element corresponding to a certain link variable on the 
rectangle, and in general 

gU(p 2 )g^ = Z{p 2 )(l N + %he%L a - l -f 2 e%e\L a L h + . . .) (3.17) 
The rectangle variable, to second order in / 1; f 2 , is 

N - ReTr[U(R)} = N(l - Re[Z( Pl )Z(p 2 )]) 

+ ±Re{Z( Pl )Z{ P2 )} (fl + f 2 + 2A/ 2 e a ■ e 3 + . . .) (3.18) 

and this is clearly as large as possible, for given fx, f 2 >0 and Re[Z \pi) Z \p 2 )] > 0, 
when ei = £3. The latter condition will always be satisfied if we choose all link 
variables in the plane to lie in the same U(l) subgroup of SU(N), so it is sufficient, 
for the purpose of proving stability, to make this choice. 
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We consider small deformations around a one- vortex configuration, which pierces 
plane at plaquette p Q , and denote the plaquettes adjacent to po, as {pi, i = 1 — 4}. 
Writing again that the deformation around a plaquette is 

V(p)=exp[if(p)e-L] (3.19) 

and defining 

_Cl 

~~ c 

£ = cos { — - ^ (3.20) 



we find the action in the plane S p i ane , up to second order in the deformation, to be 

Q EE -(Vne-iV(l-^)(l-4r)) 
Co v 

= E/ 2 (p)-^EE(/(p) + /(p + /^)) 2 

V P M=l 

+ (l-z)r^(/(p ) + /fe)) 2 - (l-z)f 2 (p ) (3.21) 

i=l 

where p + fl denotes the plaquette adjacent to p in the fi, = 1 or 2 directions tangent 
to the plane. A little rearrangement brings this to the form 

Q = £(1 - 8r)/ 2 (p) + r £ £(/(p + /2) - /(p)) 2 
p p m 

+ (1 - z)rY,(f(p ) + f{ Pl )f ■ (1 - z)/ 2 (p ) (3.22) 

i=l 

Now if z = 1, so the expansion is around the trivial vacuum, then we can 
immediately read off that Q is minimized at f(p) = 0, and the trivial vacuum is 
stable, providing that 8r < 1, which is the condition (|3.14|) above. If, on the other 
hand, 8r > 1, then Q has no minimum to second order in the deformation, and the 
trivial vacuum is unstable. 

Next, separating out of the first two sums in Q the contributions containing the 
plaquette p , we have 

Q= ]T(l-8r)/ 2 (p)+r £ (f(p + fi) - f{p)f 

+r - f( P o)) 2 + r(l - z) jjfipi) + f(p )) 2 + (z- 8r)f 2 (p ) 

i=l i=\ 

= J2(l-Sr)f 2 (p) + r (f(P + fi)-f(p)) 2 

P^PO P,P+7&P0 

+rz J2(f(Pi) ~ f(Po)) 2 + 2r(l - z) E / 2 fe) + *(1 - 8r)/ 2 (p ) (3.23) 

i=l i=l 
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If z > 0, which was the stability condition for vortices in the Wilson action, and if 
r < 1/8, which is the stability condition for the trivial vacuum in the two-parameter 
action, then every term contributing to Q is positive semi-definite. The minimum 
Q = is only obtained at f(p) = 0, i.e. at V(p) = I N everywhere, and the vortex is 
therefore a stable local minimum of the two-parameter action. 

So it appears that the result first obtained by Bachas and Dashen is quite robust: 
center vortices at N > 4 are stable minima of lattice actions in a large region of 
coupling-constant space associated with improved actions, and it seems unlikely that 
adding a few more contours, as in certain proposed perfect actions [27,28], would alter 
this result. 3 If in fact center vortices remain as local minima of the action all along 
the renormalization trajectory, then their physical effects must become apparent at 
some scale. The argument is simply that the Boltzmann suppression factor of a 
vortex in an effective (or perfect) action at lattice spacing d goes like 



cxp 



Vortex Area 



(3.24) 



where the vortex area is in lattice units. On the other hand, the entropy factor 
increases with vortex surface area as 

exp [+c x Vortex Area] (3.25) 

where c is a constant. As d increases, n 2 (d) also increases. Eventually entropy wins 
over action, and vortices at that scale will percolate through the lattice. 

So far it would appear that vortices are only stable at N > 4. At this point we 
will digress briefly to consider what happens for N = 2,3,4, following closely the 
discussion in ref. [29]. The salient point is that the effective long-range action at the 
color-screening scale should contain contours in the adjoint representation, multiplied 
by coefficients falling with a perimeter rather than an area law. The reasons for this 
are easiest to see in the context of the strong-coupling expansion. Let us start with 
the strong-coupling Wilson action with lattice spacing a, and construct a blocked 
action with lattice spacing La 

e -s'[u] = f DU Y[8[Uv - (UUU...U) l ,}e- Sw[u] (3.26) 
J v 

One readily finds the leading terms in the blocked action 
-S'[U] = £ | 2 iV (-^ ReTr F [U(p')] 

/ i \4(4L-4) 

+A(D ~ 2) [Tn) TiaW] 

+larger adjoint loops with e~ p coefficients (3.27) 



3 We have not investigated those actions systematically, however. 
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Figure 1: Diagrams leading to perimeter-law falloff, in the blocked action, of adjoint- 
representation loops. 

The adjoint loops are non-planar contributions, and arise from the "tube" diagrams 
shown in Fig. |l[ The larger adjoint loops are important for the color-screening of 
zero N-ality Wilson loops, and cannot be ignored if accurate results for such loops 
are required. 

In fact the non-local blocked action can be expressed (with some additional 
complications) as a local action containing a number of adjoint Higgs fields [29]. For 
simplicity, however, we will simply consider the above action truncated to the leading 
plaquette terms, i.e. 

S' trun M = £{c ReTr F [W(p')] + c x Tr A [U{v')]} (3.28) 

with 

c ~ exp [ — o Area(p')] , C\ ~ exp[— 4<r Perimeter (p')] (3.29) 
Consider small fluctuations around a vortex configuration 

SLnc = ^nst + ~ £ {c cos (^) Tt f [F}] + Cl Tr A [F p 2 ,]} (3.30) 

It is clear that at large blocking L, c\ is much greater than c$. Then S' trunc , at the 
scale where c\ ^> cq, is minimized at F p i = 0. This means that at this scale, vortex 
configurations are stable even for N = 2, 3, 4. 4 

4 We would like to comment at this point on an argument against vortex stability in Yang-Mills 
theory (i.e. against any definite thickness for vortices) put forward in ref. [30], and based on an 
inequality involving a certain vortex operator. We wish to point out that the argument of that 
reference could be just as well applied to Zi lattice gauge theory at strong couplings, arriving at 
the same inequality as in the Yang-Mills case, and reaching the false conclusion that Z2 vortices also 
have no definite thickness. We believe that the problematic aspect of that analysis was to identify 
a system in which vortices are restricted to a fixed set of vortex containers, associated with the 
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This result was derived for the strong coupling Wilson action. However, from 
the renormalization-group point of view, the perfect action at the color screening 
scale (which is approximately one fermi for iV = 2) is surely a strong-coupling action 
of some sort, and the same analysis should apply. There is also an interesting study 
of vortex stability in the framework of continuum QCD, that has been carried out 
recently by Diakonov and Maul [31]. 



4. Density of Vortices at Large N 

Numerical simulations in SU(2) lattice gauge theory indicate that a dilute gas ap- 
proximation for vortices is not really valid in the case of N = 2 colors [32], and a 
"vortex liquid" picture is probably more appropriate. The limitations of the dilute 
gas approximation, as we will see, are even more apparent at large N. 

Following the notation of Del Debbio and Diakonov [11], define an "/-vortex" as 
a center vortex which, when it pierces the minimal area of a Wilson loop, contributes 
a factor z l N = exp[2nil/N}. A vortex which pierces the loop area twice in opposite 
directions, and is therefore not topologically linked to the loop, contributes the trivial 
factor z l N Zx~ l = 1- The dilute gas approximation assumes that: (i) vortex piercings 
are statistically independent; (ii) there is no limitation on the number of vortices that 
can pierce a given area; (iii) it is always possible to distinguish between piercings in 
a plane due to a single /-vortex, and piercings due to, e.g, / 1-vortices. With these 
assumptions, the probability that ni /-vortices pierce a Wilson loop is given by the 
Poisson distribution 

rT ni 

P = -^-e~ n; (A 1) 

where 

n x = p(l, iV)Area(C) (4.2) 

is the average number of /-vortices piercing the minimal area of the Wilson loop, 
and p(l, N) is the number of /-vortices per unit area piercing any given plane. In the 
dilute gas approximation, the contribution of vortex piercings to an SU(N) Wilson 
loop in the /c-representation is 

N-l oo 

w(k,N)= n Z p niM klni 

1=1 m=o 

= exp[- ( T PD (A;)Area(C)] (4.3) 



derived inequality, with the actual vortex vacuum. The restriction of vortices to lie in containers 
at fixed locations would be, for the vortex vacuum, a drastic truncation of the entropy of vortices 
due to positional fluctuations. In fact, numerical experiments with both Z2 lattice gauge theory, 
and center-projected lattices in Yang-Mills theory, have shown that confinement is due to a single 
vortex configuration which percolates through the entire lattice. The vortex container of such a 
configuration has a volume comparable to the volume of the full lattice, but this fact in no way 
implies that the thickness of the vortex is the length of the lattice. 
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with the fc-string tension 



N-l 



(k) = J] p{l, N) (l - e 2 * ikl/N ) (4.4) 



1=1 



where a PD (k) and Area(C) are in physical units. The "PD" superscript is a reminder 
that this string tension is associated with the Poisson distribution 

If one accepts the assumption of a dilute vortex gas, then it is interesting to 
ask what Casimir or Sine Law scaling implies for the vortex density p(l,N). This 
question was addressed in a recent article by Del Debbio and Diakonov [11], who find 
that the scaling laws are obtained from eq. (14.4) by the vortex densities (for / <C N) 



^(4P-l) + (W) Sille LaW 

p(l,N) = Na(l) x { (4.5) 



_^ )7r 2, 2 + O (^2) Casimir Scaling 



2(jV- 

whose leading terms are proportional to N. 

A vortex density growing linearly with N is certainly pathological, but in fact this 
behavior is not what is found in center vortex theories having the fc-scaling property. 
We will first supply two examples, namely strong-coupling Z^ lattice gauge theory 
and compact QED 3 , which will serve to illustrate this point. Next, we will derive 
explicit center vortex densities leading to fc-scaling, and show that these densities are 
proportional to 1/N, rather than N. Finally, the apparent contradiction with the 
result ( |4.5| ) is explained, by showing that finite-range correlations in the center field 
strength, combined with the /c-scaling property, are inconsistent with the assumed 
Poisson distribution (|4.1|) , from which eq. ( |4.5|) is derived. 

We begin with strong-coupling lattice gauge theory, with Wilson action 



p 

_ -27m r 
-2p cos 

p 



N 



(4.6) 



Let A(C), P{C) denote the minimal area and perimeter of loop C in lattice units. For 
D = 2 dimensions, to leading order in (3, a Wilson loop in the N-ality k representation 
(with k < N/2) has a vacuum expectation value 

/ ok\MC) 

(Z k (C))= (D = 2) (4.7) 



with corresponding string tension (in lattice units) 

<7 fc = fclog(l//3)+log(fc!) (U 
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Figure 2: A cross-section of the leading strong-coupling diagram, in D > 2 dimensions, for 
the expectation value of a k = 5 Wilson loop in strong-coupling Zjy lattice gauge theory. 
The loop is perpendicular to the page, crossing it at the open and solid circles. Each 
solid line is a cross-section of a surface of plaquettes, with surfaces partially overlapping at 
the right and left-hand boundaries. In Hamiltonian formulation, the solid horizontal lines 
represent k distinct flux tubes running between a /c-quark and antiquark. 

which has /c-scaling only for k <C (3. The result is different for D > 2 dimensions. To 
leading order in f3 we find instead 

(Z k (C)) = ^A(C)+„ k P(C) ( 4 _g) 

with a string tension 

a k = ka 1 , oi=log(l//3) {D > 2) (4.10) 

which has perfect fc-scaling for k < N/2. The difference between the D = 2 and 
D > 2 cases is that for D = 2 it is necessary to expand exp((3Z(p)) to k-th order in (3 
for every piaquette in the area bounded by loop C. This brings in a factor of (l/k\). 
By contrast, in any number of dimensions greater than two, it is possible to bring 
down a single layer of plaquettes on k distinct horizontal surfaces, thereby expanding 
exp(/3Z(p)) only to first order at each piaquette on each of these surfaces. A cross- 
section of the leading strong-coupling diagram, for the case k = 5, is shown in Fig. 
||. Overlapping plaquettes at the boundaries account for the perimeter contribution 
exp(-/i fc P(C)) in eq. Q). 

What our simple example demonstrates is that /c-scaling readily coexists with 
center vortex/center monopole confinement mechanisms. Strong-coupling Z^ lattice 
gauge theory is a theory with only center vortex and (in D > 2 dimensions) center 
monopole excitations; there are no other physical degrees of freedom in the theory. 
It is also a theory in which the string tension satisfies /c-scaling perfectly, for the 
D > 2 case where center monopoles exist. 

Let us consider the vortex density p(l, N) in strong-coupling Z N lattice gauge 
theory, with lattice spacing a. The quantity p(l, N) = p(l, N)a 2 is the probability 
that an /-vortex pierces a piaquette. As f3 — > all piaquette values are equally likely, 
i.e. p — > 1/N in this limit, and in general 

p(l,N) = -^[l + 0((3)} (4.11) 
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at strong couplings. Far from growing linearly with N, this density actually falls with 
N. Here we have a complete contradiction to the behavior ( |4.5| ) derived from the 
dilute gas approximation. In the special case of D = 2 dimensions, where plaquettes 
all fluctuate independently, the density is obtained immediately from the Wilson 
action 

1 exp[2/3cos(27rZ/A0] 



P(L A ' a 2 £^exp[2/3cos(27rm/iV)] 

1 exp[2/?cos(2vr//iV)] } 



Na 2 1 + (3 2 

which illustrates the typical 1/N dependence of vortex densities in strong-coupling 
lattice gauge theory. 

The fact that strong-coupling Z N lattice gauge theory is a counter-example to 
the dilute gas result for the vortex density, eq. Q4.5| ) , is no doubt due to the fact that 
the assumptions underlying the dilute gas approximation are not satisfied by the Z N 
lattice gauge theory. In particular, in lattice theory, no more than one vortex 
can pierce a plaquette. If a Wilson loop around a plaquette has a value exp(27ril/N) 
with 1^0, this means that one I- vortex has pierced the plaquette. As a consequence, 
a planar area A(C) in lattice units can be pierced by no more than A(C) vortices. 
There is no such restriction in the dilute gas approximation, and the sum over ni in 
eq. (|4.3| ) runs to rii = oo. 

A second example of fc-scaling by center flux is provided by compact U(l) lattice 
gauge theory in D = 3 dimensions. As explained in section 2, this example is our 
paradigm for confinement by center degrees of freedom at large N. Once again, this 
is a theory in which center degrees of freedom are the only degrees of freedom in 
the theory, since the center of the gauge group is the gauge group itself. In compact 
QED 3 there is no need to invoke the strong-coupling expansion to demonstrate k- 
scaling; instead we can apply Polyakov's monopole Coulomb gas analysis to calculate 
/c-string tensions. This analysis was carried out in ref. [9], where it was found that 
fc-scaling is obtained through the existence of k independent surfaces of electric flux 
(i.e. k separate flux tubes at any given time). As we have seen, this is the same 
mechanism which accounts for A;-scaling in strong-coupling Z^ lattice gauge theory 
in D > 2 dimensions. The formation of k separate flux tubes also accounts for 
Casimir scaling, in the large N limit, in the gluon-chain model [33]. 

The example of compact QED 3 also illustrates another important point: The 
center confinement mechanism does not necessarily depend on having center flux 
of some fixed magnitude concentrated in tubes (3D) or sheets (4D). Center flux in 
compact QED 3 , passing through a given area in a plane, can spread out through 
the volume in any manner consistent with the U(l) lattice Bianchi identity; there 
is no reason for a fixed quantity of center flux to remain concentrated in a tube 
or a sheet-like region. What is essential to the center confinement mechanism is 
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that center flux passing through neighboring regions of a plane, of sufficiently large 
uncorrelated. In Z2 lattice gauge theory, where the center flux cannot 
spread out, the lack of correlation between center flux through neighboring regions 
is accomplished by percolation of vortices through the entire lattice. In compact 
QED%, where the center flux can spread out and the flux through a plaquette can be 
arbitrarily small, disordering is accomplished via a monopole plasma. In theories 
at moderate N, presumably both vortex percolation and center monopole effects are 
in play. 

The short-range correlation of center flux is sufficient, in any Zn or U(l) theory, 
to derive the area law falloff of Wilson loops. Combined with center dominance, 
short-range center flux correlation is also is sufficient to derive the area law in SU(N) 
gauge theories. The derivation for Z^ theories, including the effective P-vortex 
theory extracted (via eq. ( |2.3|) ) from SU(N) gauge theory, goes as follows: We imagine 
dividing a plane into square regions of area A, with A taken large enough so that the 
center flux piercing different regions of the plane are uncorrelated. Then consider 
a large rectangular loop C in the plane, with minimal area Area(C) ^> A, and 
subdivide this minimal area into Area(C)/^4 adjacent regions of area A bounded by 
loops {Ci,i = 1, Area(C)/*4}. In a Zn theory we have the identity 

Z k {C)=l[Z k {C i ) (4.13) 
Then, using the fact that the flux through the regions are uncorrelated, we have 

(z k (c)) = (iiz^a)} 

i 

= U(z k (c t )) 

i 

= (Z k {C l )) ATCa(c)/A (4.14) 

so that 

a(k) = --\n[(Z k (C l )}} (4.15) 

where cr(k) is the string tension in physical units. Now we define p(l,A) to be 
the probability that a region of area A is pierced by flux 2nl/N. We make no 
distinction about how that flux is divided within the region (e.g. whether the total 
flux is provided by two vortices of type 1/2, or one vortex of type /, etc.). Then 

JV-l 

(Z k (Ci)) = £ ( 4 -!6) 
1=0 



and the string tension is 



VN-l 
L 1=0 



(4.17) 
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If a Zfq lattice gauge theory has both fc-scaling and short-range correlations in 
the flux through a plane — and these conditions hold for both strong-coupling Zn 
theory and compact QED% — then we can derive the corresponding p(l,A). From 
eq. ( |4.16| ), we see that p(l,A) is the inverse discrete Fourier transform of (Z k (Ci)) 



N-l 



iV k=0 



(4.18) 



and inserting the fc-scaling behavior 

(Z\C)) = < 

we find that 



e -ka(l)A f, < | 
e -(N-k)a(l)A k > y 



(4.19) 



p(l,A) 



1 

N 



1-7 



f] + l e - 27 ri([f]+l)i/iV 1 
— + - 



7 



■$]+lJM([%]+l)l/N 



1 

N 

1 

N 



I _ ry e -2ml/N 

(1-7 2 ) (l-(-l)'7^) 



1 — 7e 



2ml /N 



2-? 

iV 



1 + 7 2 — 27 cos 
(1-7) (l + 7-2(-iy 7 [ f ] +1 cos(^)) 



for even-integer iV 



1 + 7 2 — 27 cos 



2tW 

A' 



for odd-integer A?" 

(4.20) 



where we have defined 



7 



e~°^ A and 



AT 
~2~ 



N 
2 



AT-1 



even A" 
odd A" 



(4.21) 



It can be checked that for any N, and any cr(l) > 0, 

1 , and p(l, A) > for every I 



N-l 



(4.22) 



and one can also verify that in the cr(l) — > limit, only p(0, .4.) — > 1 is non-zero (i.e. 
no string tension means no center flux). Once again we see, in eq. (|4.20|) , the overall 
factor of 1/N multiplying the flux probability distribution, as opposed to an overall 
factor of N that would have been expected from the dilute gas result (|4.5|). 



Equation ( |4.20|) is a very general result for Z^ gauge theories. It assumes only 
/c-scaling and finite-range flux correlations. The result therefore holds for any Z^ 
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theory with these properties, and these include strong-coupling Z N gauge theory 
for D > 2, compact QED 3 at any coupling, and any other theory (such as 
the complicated, non-local P- vortex theory defined in eq. ( |2.3|) ) whose long-range 
structure is described by a monopole Coulomb gas. Our result demonstrates quite 
clearly that there is no incompatibility whatever between a center flux confinement 
mechanism and fc-scaling of the string tensions. 

4.1 Why Vortex Flux Cannot Follow a Poisson Distribution at Large N 

The center flux probabilities of eq. ( j4.20|) , and the vortex density derived from the 
Poisson distribution, eq. (|4.5| ), apparently contradict one another. The former ex- 
pression is proportional to 1/N while the latter grows linearly with N, and this dis- 
crepancy needs to be explained. We will now show that the combined properties of 
fc-scaling, and the existence of finite-range correlations among center field strengths, 
are inconsistent with the Poisson distribution ( |4.1|) at large N, from which eq. ( |4.5| ) 
is derived. 

The derivation of eq. ( |4.20| ) starts from the assumption that one can subdivide 
the plane into squares of some sufficiently large area A, such that the center fluxes 
in each area are essentially uncorrelated. In every example we have suggested, where 
center flux is responsible for confinement, there is a lower bound A m in = L 2 C to the 
area A, where L c is the correlation length for the center field strength. In strong- 
coupling Z N lattice gauge theory, A m %n = a 2 , where a is the lattice spacing, so in 
this case L c = a. In compact QED3, L c is the average monopole separation. In 
center-projected QCD (the theory defined by S e ff(z) in eq. (|2.3|)), L c must be on 
the order of the characteristic length scale of the theory, i.e. L c ~ 0(Aqc D ). Below 
these length scales, all field strength fluctuations are correlated, and a treatment 
of center flux in terms of statistically independent vortex piercings is inconsistent. 
Thus, to describe the center flux in a plane in terms of statistically independent 
fluctuations, which is certainly a prerequisite for the use of the Poisson distribution, 
it is necessary to follow the procedure outlined above: Sub-divide the plane into 
squares of area A = A m im an d identify the center flux 2nl/N in a given square with 
a "piercing" of the square by a vortex of type /. 5 The vortex density is then defined 
as 

p(l,N) = S^- (4.23) 

where p(l,A) is the probability that the center flux in the square is 2nl/N. From eq. 
( |4.17| ), and the fact that the probabilities p(\,A) sum to unity, we have 

AT-1 



= -j lo£ 



p(o,^)+^p(U) 



1=1 



5 As explained above in connection with compact QED3, the term "vortex" in this context should 
not be taken to mean that the center flux of type I piercing a square region will necessarily remain 
concentrated in a tube or sheet of constant flux 2irl/N. 
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N-l 



1- P( l > A ) U-exp 



A 



log 



i=i 



2nikl 
N 



(4.24) 



where a FD {k) is the string tension in eq. ( O ) derived from the Poisson distribution. 
Eq. (|4.24j) is the correct formula for the string tension cr(k), for uncorrelated center 



flux in areas A. However, the approach based on the Poisson distribution agrees 
with this correct answer only if all the A<J PD (k) satisfy 

Acr PD {k) < 1 all ke[0,N- 1] (4.25) 

We can now see that this condition is incompatible with the fc-scaling formula 



a{k) = ka{\) (k < 



N 



at sufficiently large N, and in fact can hold only for 

2 



(4.26) 



(4.27) 



Aa(l) 

We conclude that the following conditions 

1. fc-scaling; 

2. finite correlation length L c =>- finite lower bound for A] 

3. Poisson distribution for vortex piercings; 

are incompatible at sufficiently large N. This means that the result of ref. [11], shown 
in eq. ( |4.5| ) above and indicating that fc-scaling requires a vortex density growing 
linearly with N, is not valid for realistic theories with a finite range of center flux 
correlations. 

An even simpler argument is just to note that p(0, A) ~ l/N for fixed A = A m { n 
at large N. This means that each square region is virtually certain to be pierced by 
some finite amount of center flux (as, e.g., in the case of compact QED 3 ), thereby 
violating the diluteness assumption underlying the Poisson distribution. 



5. Conclusions 

The fc-scaling of /c-string tensions at large iV appears to be a natural outcome of a 
confinement mechanism based on center degrees of freedom. Although fc-scaling, as 
the common large-N limit of Casimir and Sine Law scaling, is certainly not unique to 
the center vortex scenario, it is still of interest to see how the property emerges in that 
context. Stated briefly, the confining dynamics at asymptotic distances are described 
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at large iV by a center monopole Coulomb gas. The /c-string tension, in that system, 
is proportional to k, because there are k independent flux tubes stretching between 
the quark and antiquark. 

Thin center vortex configurations are stable classical solutions of the SU(N) 
Wilson lattice action for any iV > 4, as shown in ref. [10] and reviewed here. We have 
extended this result to a class of two-parameter lattice actions that are motivated by 
renormalization-group and tadpole- improvement considerations. It is found that for 
all couplings such that the trivial vacuum is a stable minimum of the two-parameter 
action, thin center vortices are also stable local minima of the action. We conjecture 
that for N > 4, thin center vortex configurations are stable minima of the perfect 
lattice action all along the renormalization trajectory. If that is so, then vortex 
stability at the semiclassical level is not just a lattice artifact, but is a genuine 
feature of the continuum theory. 

Finally we have shown, by explicit examples, that the /c-scaling property does not 
imply any pathological property of vortex densities in the large-N limit, of the type 
that was suggested in ref. [11]. Those pathologies derive from the use of a Poisson 
distribution for vortex densities, which we have shown here to be inconsistent with 
/c-scaling and the existence of a finite correlation length in the large-N limit. We 
have also derived the (well-behaved) vortex density distribution leading to A;-scaling 
in the large-N limit. 

The essential feature of the center confinement mechanism, aside from center 
dominance, is the finite correlation length of center flux in a plane. Independent 
fluctuations of the center flux in finite regions guarantees an area law to Wilson 
loops of non-zero N-ality, with an asymptotic string tension that depends only on 
the N-ality. The necessary randomizing of center flux at finite distances can be 
achieved via vortex percolation, as in Z 2 or Yang-Mills lattice gauge theory, or from 
the disordering due to a center monopole plasma, as in compact Z N or SU(N) gauge 
theory at large N. At moderate values of N, perhaps as low as N = 3, it is likely 
that both effects play a role. 
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